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SUMMARY 

In this study, a reliability-based structural synthesis procedure has been 
developed to tailor laminates to meet reliability-based (ply) strength require- 
ments and achieve desirable laminate responses. The main thrust of the paper 
is to demonstrate how to integrate the optimization technique In the composite 
laminate tailoring process to meet reliability design rt qu 1 rements . The ques- 
tion of reliability arises In fiber composite analysis and design because of 
the inherent scatter that is observed In the constituent (fiber and matrix) 
material properties during experimentation. Symmetric and asymmetric composite 
laminates subject to mechanical loadings are considered as application exam- 
ples. These application examples illustrate the effectiveness and ease with 
which reliability considerations can be integrated in the design optimization 
model for composite laminate tailoring. 


NOMENCLATURE 

[A] (3 x 3) composite axial stiffness submatrix 

B design variable vector; (bi , b 2 , . . b^) 

CC3 (3 x 3) composite coupling stiffness submatrix 

[D] (3 x 3) composite flexural rigidity submatrix 

<j> ply failure criterion via modified distortion energy 

E Young's modulus 

e strain vector; (e xx> eyy ( e xy ) 

F single, real valued objective function to be minimized 

G shear modulus 

M applied moment vector; (M xx> Myy ( M X y) 

N applied membrane load vector; (N xx , Nyy» N X y) 

S strengths 

ti thickness of the i^h ply 
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X 


curvature vector; ( Xxx , Xyy, x X y> 

0i ply orientation angle for the ith ply 

PI weight density of the i^ h ply 

Subscripts : 

f,m,c,l fiber, matrix, composite and ply 

k volume ratio 

L,U lower and upper bound, respectively 

T,C tension, compression 

11 direction along the fiber 

22,33 directions transverse to the fiber 


INTRODUCTION 

There is a great potential for structural tailoring of high performance, 
lightweight structures made up of fiber composites. This is because of the 
great range of desirable properties of fiber composites and their tailoring 
capacity to individual design requirements (refs. 1 and 2). Structural tailor- 
ing Is a powerful concept in that it provides the formalism to configure lami- 
nates with several significant advantages when compared to conventional Iso- 
topic material s . 

Principal among these is the ability to tailor ply properties (material 
selection) and fiber orientation (configuration selection) to the given mechan- 
ical, thermal and hygral loads. Thus, a laminate is a building block that con- 
tributes to the overall thermostructural action of a component made out of 
composites. Hence, it is only logical to concentrate on tailoring of laminate 
properties and configuration as a first step in the synthesis of composite 
structures in order to meet reliability design requirements. 

The idea of applying optimization techniques to automate the structural 
tailoring of composite components such as engine blades with the aid of comput- 
ers is not new. The structural tailoring of engine blades (STAEBL) and struc- 
tural tailoring of advanced turboprops (STAT) computer programs have been 
developed and routinely used at NASA Lewis Research Center (refs. 3 and 4). 
Almost all of the existing structural tailoring software is based on a deter- 
ministic approach. Thus it is not possible to directly account for the inher- 
ent composites properties scatter, in the currently available tailoring 
software. To be realistic, uncertainties in the composite material properties 
must be reflected in the analysis and design methods for structural and mechan- 
ical systems. Chami s (ref. 5) has underscored the importance and usefulness of 
probabilistic approach in the analysis of space propulsion systems to Improve 
the reliability of engine components. 
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In general, the Input data prone to uncertainties during structural analy- 
sis Include material properties, boundary conditions, loading conditions and 
geometry of the system. In the present research, we shall focus our attention 
to uncertainties In the material properties and how to incorporate them In the 
laminate tailoring process. The most common and traditional approach to 
account for uncertainties in the design process is to -Introduce a factor of 
safety. But this approach does not quantify the reliability of the system 
design. Recently, probabilistic concepts in the form of Monte Carlo simulation 
have been applied to quantify the uncertainties in composite micromechanics by 
Chamis and coworkers (ref. 6). Such probabilistic composite micromechanics 
concepts suggest a viable approach towards integrating reliability considera- 
tions in laminate tailoring in the form of behavioral constraints and/or proba- 
bilistic load conditions as a first step. 

The main objective of this research effort is to develop a tailoring 
(synthesis/optimization) procedure by combining a robust optimizer with compos- 
ite mechanics and probabi ’ 1 sti c constraints and/or loads. It is widely recog- 
nized that at the optimum design, one or more performance constraints are at 
their allowable values and any uncertainties in the given input data may for 
these constraints lead to an unsafe design. This is because the active per- 
formance constraints may become violated with variations in the input data. 

The procedure described herein, is applied to several generic sample cases to 
illustrate that fiber composite laminates can be tailored to meet design 
requirements with assured probability of survival. 


PROBABILISTIC COMPOSITE MICRCMECHANICS 

The branch of composite mechanics that predi :ts ply material properties 
based on the constitutent properties, volume fraction end fiber orientation is 
known as composite micromechanics, and frequently Incorporates the traditional 
mechanics of materials assumptions. As mentioned in the introduction, Chamis 
and coworkers (ref. 6) have developed a probabilistic approach to composite 
micromechanics using Monte Carlo simulation. Thus, scatter in the primitive 
variables namely constituent (fiber and matrix) elastic properties, strength:, 
and fiber volume ratio have been quantified at the ply level properties. A 
straightforward way of developing probabilistic analysis is to perform Monte 
Carlo simulation using deterministic composite micromechanics equations. 

Thus, fiber volume ratio, fiber and matrix elastic properties and strengths 
are considered to be independent random variables. The scatter in these inde- 
pendent random variables is quantified by defining the type of probability dis- 
tribution function (PDF) e.g., normal, Weibull etc. and the related statistical 
quantities such as the mean and the coefficient of variation etc. The avail- 
able experimental data usually provide the estimation of the mean value and 
the coefficient of variation. 

The ply level properties such as moduli (longitudinal, transverse and 
shear), strengths (longitudinal tensile/compressive. Transverse tensile/ 
compressive and intralaminar shear) are considered to be dependent random var- 
iables. The results of the probabilistic analysis are presented in terms of 
(1) histogram or PDF, and (2) cumulative distribution function (CDF) of the 
ply level properties. Typical PDF and CDF curves for the longitudinal tensile 
strength (Sqjij) of a carbon graphite fiber (AS)/intermediate modulus high 
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strength (IMHS) epoxy matrix composite system ply (ref. 6) are shown in figure 
1. Note that one can readily assess the probability of a random sample to be 
higher than a given value for Sgju from the CDF curve. In this study, the 
longitudinal tensile/compressive, transverse tensile/compressive and intralami- 
nar shear strengths of an AS/ IMHS ply are assumed to be random variables and 
their probability distribution function curves for case 2 from reference 6 
have been utilized. Table I gives the deterministic mean values of the 
strengths (Sqjit, S&11C* S&22T* Sq. 22C> ^8.12S^ a unidirectional AS/ IMHS ply 
along with the strengths corresponding to 90 percent cumulative probability. 


STRUCTURAL TAILORING SOFTWARE 

The computer program IDESIGN (Interactive Design Optimization of Engineer- 
ing Systems) is an interactive general purpose optimization software suitable 
for tailoring of laminate properties (ref. 7). The IDESIGN program has been 
widely used in static structural, dynamic, distributed parameter, shape optimi- 
zation and numerous other application problems (refs. 8 and 9). In summary, 
IDESIGN program solves the following nonlinear programming problem (NLP) with an 
objective function (F) to be minimized subject to the equality and inequality 
constraints. Thus, the design optimization model, in the form of an NLP prob- 
lem, Is to find a k-dimensional design variable vector, B to: 

Minimize F(B) such that 

G(B) <0 (1) 

B L < B < Bu 

Note that in the above NLP problem defined by equation (1), it is assumed that 
the functions F and G are continuous and differentiable. The design varia- 
bles (B) are also continuous real variables. The Sequential Quadratic Program- 
ming (SQP) method is chosen as the optimization algorithm (refs. 10 and 12). 

The numerical data such as starting values and upper/lower bounds on the design 
variables, number of design variables and constraint , tolerances on constraint 
violation and convergence criterion are provided to .he IDESIGN program through 
an input data file. The gradients of the objective and constraint functions 
are calculated by a forward finite difference scheme in the IDESIGN program. 

In every design iteration of the SQP method, the following quadratic program- 
ming (QP) subproblem is solved to get the search direction d at a design 
point B, 


Minimize { F ( B) + VFTd + 0.5 d T Hd) 
such that G(B) + VG^.d < 0 


( 2 ) 


where VF and VG are gradient vectors of the objective and constraint func- 
tions; Bp is the approximation to the Hessian of the Lagrange function. The 
details about the QP subproblem and its solution can be found in Ref. 9. A 
suitable step size (a) needs to be found along the search direction d so 
that progress can be made towards the optimum solution. Note that the search 
direction (d) is a vector whereas the step size (a) is a scalar. Thus, the 
change in design is given as a • d and the recursive formula for the next 
design point i s given as, 
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( 3 ) 


The SQP method converges when the residue, ||B new - B°^|| becomes very 
small. Overall, IDESIGN software is interactive, user friendly and modular. 

Any application/analysis program can be easily coupled with the IDESIGN system. 

The X n tegrated Composite Analyzer (ICAN), a stand-alone computer code, has 
been used to analyze and design multilayered fiber composite laminates using 
micromechanics equations and laminate theory (ref. 13). Input parameters of 
this user-friendly program include material system, fiber volume ratio, lami- 
nate configuration, fabrication factors, and environmental conditions. Output 
features Include the practically most of the composite hygral , thermal and 
mechanical properties that are required to perform structural /stress analyses 
in service environments. As such, ICAN may be utilized as an effective tool 
for the preliminary design of composite structures (ref. 14) as shown in fig- 
ure 2. In addition, ICAN has a resident data bank which houses the properties 
of a variety of constituent (fiber and matrix) materials with provisions to add 
new constituent materials as they become available. Note that ICAN program Is 
converted into a subroutine for the IDESIGN software to perform synthesis of 
laminate properties. Thus, the ICAN code is utilized to evaluate the objective 
and constraint functions in equation (1). The overall flowchart for the compu- 
tations performed In the proposed structural tailoring program is given in 
figure 3. 


LAMINATE CONFIGURATION 

An eight-ply laminate configuration made up of AS-IMHS composite material 
system Is considered herein. Both symmetric and asymmetric laminate configura- 
tions are considered. The use and cure temperatures are both 70 °F with 
0 percent moisture content. The individual ply orientation angle, ply thick- 
ness and fiber volume ratio for various plies in a laminate can be defined in 
the ICAN input data file. 


PROBLEM DESCRIPTION 

The laminate tailoring tasks considered herein can be grouped under two 
categories: 

Category 1: The application examples in this category are deterministic in 

nature and are mainly selected to test the laminate tailoring approach proposed 
in the paper. Here, no hygral, thermal or mechanical loads are applied on the 
laminate. The eight-ply laminate is confined to be symmetric. The design var- 
iables are constituent (fiber and matrix) material properties, fiber volume 
ratio (kf) and the ply angles (0^). Thus, the total number of design variables 
is 34. Only one composite material system AS-fiber in an intermediate-modulus- 
high-strength matrix (AS/IMHS) is considered in this example. There are no 
performance constraints. 

The following laminate tailoring application examples were studied in 
category 1 : 
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( 1 ) To simultaneously minimize coefficients of thermal expansion (a) and 
maximize thermal conductivities (b) and maximize shear modulus (c). Note that 
a four-ply laminate Is considered for this example. The classical optimization 
model cannot deal with multi-objective function problem. Thus, one needs to 
construct a composite objective function so as to minimize (a + 1 /b + 1/c), 
where a, b, and c denote coefficient of thermal expansion, thermal conduc- 
tivity, and shear modulus, respectively. 

( 2 ) To simultaneously maximize the elastic moduli along laminate struc- 
tural axes x, y, and z. Thus, the composite objective function is to minimize 

E cxx + E cyy G cxy 

E Slll E 2.22 G (ll 2 

where (E^ii), C Eq, 22 ^ an d ^ ^0.1 2 ^ are the corresponding values of moduli of a 
unidirectional laminate in the material axe, (1, 2, and 3). 

( 3 ) To maximize the twisting stiffness (D33) while constraining the cou- 
pling between bending and twisting represented by 0^3 and D23 terms in the 
flexural rigidity matrix (D) to zero. Note that the force displacement rela- 
tionships for a composite laminate are given as: 

-n 

Thus, the problem is to maximize D33 such that D13 = 0 and D23 = 0 . 

( 4 ) To minimize the coupling between bending and twisting actions by mini- 
mizing D] 3 and D23 terms of the flexural rigidity matrix in equation ( 5 ) 
simultaneously. This shall minimize the asymmetry in the laminate configura- 
tion. The objective function is to minimize < 3 + D23). 

( 5 ) To maximize the bending as well as twisting stiffness of the lami- 

nates, or minimize the negative of the diagonal elements in the flexural rigid- 
ity matrix i.e., to minimize: (- D] 1 - D22 - D33). 

(6) To maximize the extensional stiffness of the laminate, or minimize 

the negative of the diagonal elements in the axial stiffness matrix (A) in 
equation ( 5 ) I.e., minimize: (- A]] - A22 - A33). 

( 7 ) To minimize the shear stretch coupling effect in a laminate, or mini- 
mize the off diagonal elements of the axial stiffness matrix I.e., minimize: 

< A 1 3 + A 2 3 >. 

Category 2 : The application examples in this category are probabilistic 

In nature. They represent point designs for specified cumulative probability 
of individual strenghts. Here, mechanical loads (axial forces and bending 
moments) are applied to the eight-ply composite laminate. The design variables 
Include fiber volume ratios, ply orientation angles and ply thicknesses. Thus, 
the asymmetric laminate has 24 design variables whereas the symmetric laminate 
has 12 design variables. The constraints are Imposed on :he ply failure crite- 
ria such that the modified distortion energy (for combined stress failure) is 


r°> 

) X l 
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always positive for each ply. Thus, there are eight constraints. The modi- 
fied distortion energy <<J>) is defined in the ICAN program as: 




\2 

q Q,223 ] _ k CT ftlla 

S Q,22|V 2,12 S !llla S S122|3 



( 6 ) 


The ply failure criterion constraint is written as, 

<j> > 0 

The reader is referred to ICAN user's and programmer's manual for the details 
of equation (6). Note that the ply strengths S in equation (6) to be used 
are the values corresponding to 90 percent cumulative probability as given in 
table I. The corresponding probabilistic ply constitutent properties and any 
small deviations in the ply orientation angle are included in the scatter 
assumed for the fiber volume ratio, 0.45 to 0.55, which was used to develop the 
CDF curves for strength in figure 1. The objective function to be considered 
Is to minimize the weight per unit area of the composite laminate i.e., 

8 

E oih 

i.i 

where pi and t^ are the weight density and the thickness of 1^ ply, 
respectively. 

The following upper and lower bounds are Imposed on the fiber volume 
ratios, ply angles and thicknesses for all the plies, 

-90° < 0i < 90° J 

0.45 < (k f ). 0.55 > 1-1,8 

0.005 in. < t^ < 0.55 in.j 

The initial ply stacking sequence for all the bases is (2S/-25/25/-25/— 25/25/ 
-25/25), the fiber volume ratio and ply thickness is 0.45 and 0.01 respec- 
tively, for all the plies. The following example problems are studied under 
the probabilistic category: 

(8) To minimize the weight of the eight-ply composite asymmetric laminate 
subject to the satisfaction of the ply failure criterion for all the plies 

= 1000 lb/in., are applied, 
al constraints. 

(9) This example is the same as example 8 except that the laminate is 
considered to be symmetric. Thus, there are 12 design variables and 8 
constrai nts . 


where axial loads, N xx =* 3000, Ny y = 1000, and N xy 
This example has 24 design variables and 8 behavior 
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(10) To minimize the weight of the eight-ply symmetric composite laminate 
subject to the satisfaction of the ply failure criterion for all the plies 
where bending moments, M xx = 100, M yy = 100, and M xy = 100 lb in. /in., are 
appl ied. 

(11) To minimize the weight of the eight-ply symmetric composite laminate 
subject to the satisfaction of the ply failure criterion for all the plies 
where axial loads, N xx = 3000, N yy = 1000, N xy = 1000 lb/1n., and bending 
moments, M xx = 100, M yy = 100, M xy = 100 lb in. /in., are applied. 


RESULTS AND DISCUSSION 

The results for the example problems in the first category (determinis- 
tic) are summarized In table II. As expected, the ply orientation angles and 
fiber volume ratio are the dominant design variables in the tailoring process. 
Note that less than 25 design cycles are required for these example problems 
to converge to the optimum designs. The final results obtained in the example 
problems agree with the intuitive reasoning guided by the physics of the prob- 
lem and appear to be superior than what can be achieved manually by using para- 
metric studies. 

The ply orientation angle, fiber volume ratio and ply thickness for earh 
ply for the examples in the second category are given in tables III to VI for 
example numbers 8 to 11, respectively. The initial and final values for the 
weight of the composite laminate is given in table VII. The ply stress con- 
straints for the Initial and final values are summarized in table VIII. Note 
that the starting design for all the example problems was infeasible because 
the initial design violates the probabilistic ply strength constraints. 

The point probabilistic designs described In cases 8 to 11 can readily be 
extended to generate cumulative distribution functions for probability of fail- 
ure versus weight. Three such examples are described below. 

(1) Probabilistic distribution of the ply strengths (S^nt* Sqji c , S&22t> 
S&i2s) while loads are assumed to be deterministic. 

(2) Probabilistic distribution of the loads on the laminate whereas the 
ply strengths are fixed to a predetermined level of probability of failure. 

(3) Combined probabilistic load and strength. 

In case (1), the ply strengths are selected for differen probability lev- 
els of 50, 60, 70, 80, and 90 percent from the respective CDF curves as shown 
in figure 1 for longitudinal tensile strength (S^nt). The laminate weight per 
unit area is then optimized subject to each of these probabilistic strength 
constraints for the above mentioned probability levels. The optimum laminate 
weight corresponding to different strength r e 1 i ab i 1 ' ty levels can then be plot- 
ted (fig. 4). Note that the laminate loads are assumed to be deterministic 
and thus fixed. 

In case (2), the laminate loads are randomly selected for different proba- 
bility levels between the range 80 to 120 percent from the a probability 
density curve. The mean values for the loads are N cxx = 3000, N cyy = 1000 and 
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N CX y = 1000 lb In. /In. Once again, the laminate weight Is optimized subject 
to the strength constraints for a given probability and subject to the probabi- 
listically varying loads. The optimum laminate weight corresponding to differ- 
ent load probability levels can then be plotted (fig. 4). 

In case (3), both load and strengths are probabilistically described. The 
resulting curve is shown in figure 4. These three cases, collectively, demon- 
strate how optimization techniques can be routinely used to tailor laminates 
for probabilistic design requirements. The nearly flat curve for probabilistic 
load with fixed probability for strengths lends credence to point probability 
design. Stated differently, the optimizer will select the design variables to 
meet a specified reliability. For example, in a specific design, the struc- 
tural reliability is assured after the probability for strength excedence has 
been selected. This is independent of the loading conditions and their 
respective scatter. 


CONCLUSIONS AND SUMMARY 

A procedure has been developed for the tailoring of composite laminates 
subject to probabilistic behavior constraints and/or loads. The probabilistic 
nature of the composite material properties such as ply strengths is consid- 
ered by constraining the properties to a known probability level on the cumula- 
tive distribution function (CDF) curve. The CDF curves for the ply mechanical 
properties are obtained from probabilistic composite micromechanics. One of 
the main advantages of the procedure is that the probability of failure need 
not be evaluated in the laminate tailoring process. The example prob'ems 
described demonstrate how to tailor composite laminates by accounting for 
uncertainties at the material and/or loads level. It is concluded from the 
example cases that laminate tailoring with probabilistic constraints and/or 
loads for point and/or distributed designs can be integrated in existing struc- 
tural tailoring software for composites in a simple and straightforward manner. 
The results show that simultaneous consideration of probabilistic loads and 
strength constraints results in optimum designs with lowest reliability. 
Probabilistic strength constraints yield optimum designs with the highest 
rel lability. 
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TABLE I. - PLY VALUES OBTAINED FROM 
DETERMINISTIC AND PROBABILISTIC 
ANALYSIS 


Ply property 

Determini sti c , 
ksi 

Probabi 1 i sti c a , 
ksi 

E £l 1 

15 750 

— 

E <?22 

1 065 

— 

G«12 

516 

— 

$01 1 T 

203 

150 

$<21 1C 

165 

85 

S 022T 

11.74 

4 

S 022C 

27.41 

9 

S«12S 

10.01 

5.5 


d The property values correspond to 90 percent 
cumulative probability on the cdf curves 
given in reference 6 for case 2. 


TABLE II. - RESULTS FOR EXAMPLES IN CATEGORY I 


Example 

Design objective 

Initial design 

Final design 

1 

Minimize: coefficient of 
thermal expansion 
Maxmimize: thermal 
conduct i vi ty 
Maximize: shear modulus 

( 25/-2S) s 
k f = 0.45 

( 30 . 9/-30 . 9) s 
k f = 0.65 

2 

Maximize: E c ^ 
E cl 2 

e c22- 

(25/-25/25/-25) s 
k f = 0.45 

(-13/13/-1 3/13)- 
k f = 0.65 

3 

Maximize: D 33 
°13 : 

- 0 , D 2 3 - 0 

(25/-25/2S/-25) - 
k f = 0.45 

(0/-83/68/-33) - 
k f = 0.^5 

4 

Minimize: D^, 

°23 

(25/-25/25/-2S) s 
k f = 0.45 

( 0 / 0 / 0 / 0 ) s 
k f = 0.25 

5 

Maximi ze: D^ ^ , 

d 22 . D 33 

(25/ —25/25/— 25 ) s 
k f = 0.45 

( 0 / 0 / 0 / 0 ) s 
k f = 0.65 

6 

Maxi mi ze: &]] , 

a 22 , a 33 

(25/-2S/25/-25) s 
k f = 0.45 

( 0 / 0 / 0 / 0 ) s 
k f = 0.65 

7 

Minimize: A^, 

A23 

(25/-2S/25/-25)- 
k f = 0.45 

(-90/90/-90/90)- 
k f = 0.25 
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TABLE III. - RESULTS FOR EXAMPLE 8 a 


8 

Problem: minimize ^ p.t. such that • 4>(I) > 0.01; I = 1,8 

i = 1 

N = 3000, N = 1000, N = 1000 
xx yy xy 


Ply 

1 

2 

3 

4 

5 

6 

7 

8 

9. 

i 

26 

-26 

35 

-35 

-35 

35 

-26 

-26 

t. 

l 

0.005 

0.005 

0.005 

0.005 

0.045 

0.055 

0.0055 

0.0094 

<V. 

l 

.55 

.53 

.45 

.45 

.45 

.45 

.45 

.55 


a Unsymmetric ply angle layup. 


TABLE IV. - RESULTS FOR EXAMPLE 9 
8 

Problem: minimize p.t. such that • 4>(I) > 0.01; I = 1,8 

ft 1 1 

N = 3000, N = 1000, N = 1000 
xx yy xy 


Ply 

1 

2 

3 

4 

5 

6 

7 

8 

0 . 

i 

27 

-27 

35 

-35 

-35 

35 

-27 

27 

t. 

i 

0,0068 

0.0067 

0.005 

0.005 

0.005 

0.005 

0.0067 

0.0068 

(k f ) 

i 

.53 

.53 

.45 

.45 

.45 

.45 

.53 

.53 
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TABLE V. - RESULTS FOR EXAMPLE 10 


8 

Problem: minimize ^ p.t. such that • <|>(I) > 0.01; I = 1,8 

i = l 1 1 

M = 100, M = 100, M = 100 
xx yy xy 


Ply 

1 

2 

3 

4 

5 

6 

7 

8 

9. 

l 

45 

-25 

25 

-25 

-25 

25 

-25 

45 

t. 

i 

0.023 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.023 

<V. 

i 

.45 

.45 

.45 

.45 

.45 

,45 

.45 

.45 


TABLE VI. - RESULTS FOR EXAMPLE 11 


8 

Problem: minimize ^ P.t. such that * 4>(I) > 0.01; I = 1,8 

i = l 

N = 3000, N = 1000, N = 1000; M = 100, M 100, M = 100 
xx yy xy xx yy xy 


Ply 

1 

2 

3 

4 

5 
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00 

0. 

i 

45 

-25 

25 

-25 

-25 

25 

-25 

45 

t. 

i 

0.027 

0.017 

0.005 

0.005 

0.005 

0.005 

0.017 

0.027 

( V. 

i 

.45 

.55 

.45 

.55 

.55 

.45 

.55 

.45 


TABLE VII. - INITIAL AND FINAL 
LAMINATE WEIGHT WITH 
PROBABILISTIC 
CONSTRAINTS 


Examp! e 


Optimum 


B H 






.0042 

.0051 

11 

.0042 

.0058 
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TABLE VIII. - LAMINATE TAILORING WITH PROBABILISTIC CONSTRAINTS INITIAL (I) AND FINAL (F) VALUES 

(a) Example 8 
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flGURE 5. - Ft FWCHART TOR COMPOSITE I AMINAFE TAILORING WITH PROUARI i I STIC CONSTRAINTS. 
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